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We study relativistic charmonium on very fine quenched lattices (/3 = 6.4 and 6.6). We concentrate on the 
calculation of the hyperfine splitting between r/c and J/iii, aiming for a controlled continuum extrapolation of this 
quantity. Results for the r?c and J /ill wave functions are also presented. 



1. Introduction 

Charmonium spectroscopy on the lattice is not 
at all straightforward. Charm is too heavy for 
most current simulations - typically ma > 1 ~ but 
too light to blindly rely on the heavy quark ap- 
proximation - for cc, I (? ~ 0.3. A good probe 
of relativistic effects is the hyperfine splitting be- 
tween the ■^S'l and the ^5*0 states, which for char- 
monium is AM = 117 MeV. Lattice quenched cal- 
culations underestimate AA/ by about 30 — 50% 
(l|||l, the prediction from NRQCD being AAf = 
55(5) MeV The discrepancy could be due 
to the quenched approximation. However, first 
estimates including dynamical quarks seem to in- 
dicate that they account for only ~ 10% of the 
difference Q. Here we compute AAf on fine lat- 
tices within the relativistic formalism, aiming for 
the extraction of a controlled quenched contin- 
uum limit. 

Two approximations will be made: (i) 
quenched, (ii) OZI, meaning that Zweig-rule for- 
bidden diagrams, although contributing to singlet 
mesons like charmonium, will not be included. 



*Talk presented by M. Garcia Perez 



2. Results 

Configurations are generated with the standard 
Wilson action on a 32'^x96 lattice at /3 = 6.4 and 
6.6. Quark propagators are computed using Wil- 
son, tree-level clover and tadpole improved clover 
Dirac operators (at = 6.6, c^t'^P = 1.388 while 
the non-perturbative Csw — 1.467 1^). We have 
60 configurations at each /3 value. 

2.1. Hyperfine splitting 

AA/ is computed from the ratio of vector to 
pseudoscalar correlators. Fig. |l| shows the ef- 
fective mass at /3 = 6.6 for a pseudoscalar mass 
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Figure 1. Hyperfine splitting effective mass. 
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tuned close to the physical rjc mass (the scale set 
by ro from [^). Even on this fine lattice, the sen- 
sitivity of AM to the Dirac operator is large, tad- 
pole improvement giving the best estimate. The 
total set of results for the pseudoscalar mass and 
the hyperfine splitting is presented in Table |l|. 

Table 1 

Pseudoscalar mass (Mps) and hyperfine splitting 
(AM) in lattice units. W, CI and Tadp denote 
Wilson, clover and tadpole improved clover oper- 



ators. 


Op. 




K 


aMps 


aAM 


w 


6.4 


0.1389 


0.7632(8) 


0.01008(30) 


CI 


6.4 


0.1324 


0.7645(8) 


0.01664(44) 


w 


6.6 


0.1415 


0.5228(9) 


0.01176(23) 


w 


6.6 


0.1400 


0.5993(9) 


0.00915(24) 


w 


6.6 


0.1385 


0.6729(9) 


0.00763(19) 


w 


6.6 


0.1375 


0.7205(9) 


0.00665(17) 


CI 


6.6 


0.13225 


0.6735(8) 


0.01157(35) 


CI 


6.6 


0.1335 


0.5969(9) 


0.01357(34) 


Tadp 


6.6 


0.1158 


0.6700(9) 


0.01371(34) 


Tadp 


6.6 


0.1167 


0.6031(9) 


0.01538(38) 



To perform the continuum extrapolation we use 
13 = 6.0 and 6.2 data from UKQCD 1,0], available 
for both Wilson and tree-level clover Dirac opera- 
tors. Fitting the dependence of AM versus 1 /Mps 
at fixed f3 allows to interpolate AM(/3) to a com- 
mon value of the pseudoscalar mass: Mps = 2.93 
GeV, which is our value at /3 = 6.4 (experimen- 
tally M^^ =2.98 GeV). In Figures |, | we present 
the continuum extrapolation of the Wilson and 
tree-level clover hyperfine splittings, respectively. 
Note that: («) A linear extrapolation in a from 
f3 = 6.0 is not justified. This is a clear indication 
of strong lattice artifacts on this coarse lattice 
(aMps ~ 1.5). (m) a continuum extrapolation is 
mandatory: at all the simulated /? values the dif- 
ference between Wilson and clover data is large 
but their extrapolations are consistent with each 
other. We quote, from the quadratic (linear) fits, 
AM = 87±4 (75±4) MeV and AM = 99±7 (87±2) 
MeV for Wilson and clover respectively. We have 
also included in Fig. ^ the value of AM obtained 
from tadpole improvement at /3 = 6.6. A more 
extensive calculation including non-pcrturbative 
improved clover data is in progress. 
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Figure 2. Continuum extrapolation of AM with 
the Wilson Dirac operator. 
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Figure 3. Continuum extrapolation of AM with 
the tree-level clover improved Dirac operator. 

2.2. Wave functions 

Within the naive non-relativistic approxima- 
tion it is easy to understand the origin of the 
hyperfine splitting. This approximation amounts 
to solving the Schrodinger equation in a non- 
relativistic potential and dealing with relativis- 
tic corrections in perturbation theory. To zeroth 
order '^Si and ^Sq states are degenerate. The 
first order correction to the mass is the aver- 
age of the spin-spin interaction giving AM oc 
Q;s(mg)|^'(0)p/m^, with vE'(O) the value of the 
non-relativistic wave function at the origin. Per- 
turbative corrections to the wave function also 
depend on the spin; to lowest order, |vE'(0)| in- 
creases (decreases) for the pseudoscalar (vector). 

We have extracted gauge invariant wave func- 
tions from lattice 4-point functions: 5'(.t) = 
{cc\{^l^c){0){'4's'tpc){x)\cc) . A comparison be- 
tween pseudoscalar and vector wave functions is 
presented in Fig. ^. The observed pattern corrob- 
orates qualitatively the predictions of the heavy- 
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Table 2 



Hyperfine splitting from different approaches. The scale is set by rp or the spin averaged ^Pi—lS splitting. 





This work 


This work 


CP-PACS D 


CP-PACS |1|1 


latest NRQCDQ 




Relativistic 


Relativistic 


Relativistic 


Relativistic 


Non relativistic 




Isotropic 


Isotropic 


Anisotropic 


Anisotropic 






Continuum 


Tadpole (3 = 6.6 


Continuum 


Continuum 




Scale 


''0 


ro 




iPi - IS 


iPi - IS 


AM(MeV) 


83-106 


78(2) 


65(1) 


82(2) 


55(5) 



quark model. We also exhibit the difference be- 
tween the wave functions extracted from Wil- 
son, clover and tadpole improved clover Dirac op- 
erators. As already observed for the hyperfine 
splitting, relativistic effects (implying degeneracy 
breaking between ^Si and -^Sq) increase as we im- 
prove the fermionic action. 
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Figure 4. Pseudoscalar (top) and vector (bottom) 
wave functions. 

2.3. Conclusions 

Table I compiles data from We compare 

with an interval covering our Wilson and clover 
continuum extrapolations. Using in all cases tq 
to fix the scale, our relativistic extrapolation gets 
much closer to the experimental value. Note that 
this is also the case for our f3 — 6.6 tadpole im- 
proved clover data, which gives AA/ = 78(2) MeV. 



We conclude that a careful and controlled 
quenched continuum extrapolation is now feasi- 
ble, and gets closer to the experimental value 
than previous determinations. Dynamical quark 
effects, which presumably account for the bulk of 
the remaining disagreement with experiment, ap- 
pear to be 0(10-20)%, which is consistent with 
the change in the coupling constant a(A^/ = 0)— > 
q;(A^/ = 2) entering in the non-relativistic expres- 
sion. OZI suppressed contributions are expected 
to be small for these heavy quark masses, but we 
cannot rule out the possibility that they give a 
non negligible contribution to AM, which could 
be enhanced if there is mixing with glueballs. 
Work is in progress to improve the extrapolations 
and to extract masses for other channels. 
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